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ORBIT TRIM PROPULSION REQUIREMENTS FOR 
SUN SYNCHRONOUS SATELLITES 


Daniel L. Endres 


ABSTRACT 

Two problems are considered in this report. The first is to 
determine what method of thrusting would be the best for correct- 
ing the precession rate of a near sun synchronous satellite to the 
precession rate of a sun synchronous satellite. The second is to 
compare five propulsion systems that could fulfill the requirements 
of such a mission. Special emphasis is placed on the Tiros M 
mission. 

It is found that a tangential thrust is the best means for cor- 
recting to a sun synci ronous orbit from a propulsion standpoint. 
The results of the propulsion system comparisons are that: 1) the 
2 lb. hydrazine system and the 200 sec., I sp resistojet weigh about 
the same, 2) the 2 lb. hydrazine system requires the least time and, 
3) the cold gas resistojet and the 2 lb. hydrazine system require 
very little power. Thurst misalignment is not considered. 
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LIST OF SYMBOLS AND CONSTANTS 
I = Inclination, degrees 

0 = Average precession rate, radians/second unless otherwise specified 

J 2 o - Second zonal harmonic in the series expansion of the gravity potential 
for the Earth, 1.082 x 10" 3 

T) = Orbital rate, radians/second 

r e = Mean equatorial radius of the Earth, 3440 nautical miles 

r = Satellite radial distance measured from the center of the Earth, ft. 
unless otherwise specified 

/x = Gravitational constant of the Earth, 62,746.8 n.mi. 3 /sec. 2 or 
1.4076449 x 10 16 ft. 3 /sec. 2 

M = Moment vector, ft.-lbs. 

H = Angular momentum, slug.-ft. 2 /sec. 

V = Angle, in the orbit plane, between the line of nodes and radius vector, 
radians 

/S A A 

i , j . k = Unit vectors along X, Y, and Z respectively 

M p = Precession moment due to oblateness of earth, ft.-lbs. 

M t = Moment due to thrust, ft.-lbs. 

F n = Normal thrust, lbs. 

F t = Tangential thrust, lbs. 
m = Mass of satellite, slugs 
T = Period of orbit, sec. 

I T = Total impulse, lb.-sec. 

W 0 = Initial satellite weight, lbs. 


xi 



x -p 

av t 

AV„ 


'e f f 


Ar 

R i 

r 2 

g 

T T 

r, 6, <p 
Sr, S$, S0 

F r . F # . F. 
n 
r 


I 


mean 


L 


a; 


7 


Specific impulse, sec. 

Total effective velocity increment for tangential thrust method, 
ft./sec. 

Total velocity increment for normal thrust method, ft./sec. 
Av n /av t = penalty factor 

Final sun synchronous altitude minus reference altitude, n.mi. 
Radius of initial circular orbit, ft. 

Radius of final circular orbit, ft. 

Acceleration due to gravity, ft./sec. 2 
Thrust on time, sec. 

Polar coordinates 

Perturbations of polar coordinates 

Forces in r, d t and <p directions respectively, lbs. 

Constant angular velocity, rad./sec. 

Dimensionless time 

Mean change in inclination, radians unless otherwise specified 
Change in inclination over an orbit, radians 
Longitude of satellite with respect to ascending node 
Latitude of satellite 

Angle between ascending node and perigee 
True anomaly of satellite, radians 
Satellite orbital rate, rad./sec. 

Oblateness force acting on satellite, lbs. 


xii 


Ah 

'p. 

W o 

w t 


W 


th 


w 


f s 


W PC 

W tel 

Kj ( i=0, 3) 

cr v a 2 


Change in altitude, n.mi. 

Weight of propulsion system, lbs. 

Propellant weight, lbs. 

Tankage weight, lbs. 

Weight of thruster, lbs. 

Weight of feed system, lbs. 

Weight of power conditioner, lbs. 

Weight of telemetry, lbs. 

Structural constants for the propulsion systems 
Conversion factors 


xlii 



BLANK PAGE 


ORBIT TRIM PROPULSION REQUIREMENTS FOR 


SUN SYNCHRONOUS SATELLITES 


INTRODUCTION 

There are two main purposes to this report. The first is to determine the 
thrust requirements necessary to adjust an initial near sun synchronous orbit 
to a final orbit which is sun synchronous. It will be assumed that the satellite 
is in a near-polar and near-circular orbit about a polar oblate Earth. The 
second is to compare the weights, the orbit correction times, and the power re- 
quirements of various propulsion systems that could fulfill the thrust requirements. 

An obvious advantage of a near polar sun synchronous orbit is that the power 
output can be optimized. This is due to the larger amounts of radiant energy that 
can be encountered by the solar panels. 

In various reports and texts it has been shown that a near circular satellite 
orbit subjected to the polar oblateness of the earth’s gravity field will precess 
at a rate that depends on the orbit's semi-major axes and inclination (c.f. 
References 1 and 2). This then leads to the possibility of defining a reference 
orbit in which the satellite will precess at a sun synchronous rate without the use 
of any thrust when not subjected to other external orbit perturbations. Such a 
reference orbit can in fact be defined by an equation which is derived in Appen- 
dix A. 

If the satellite is not put into a reference orbit, thrust can be applied to ob- 
tain the proper precession rate. There are two general ways this can be done: 

1) use a continuous thrust to increase or decrease the precession that is due to 
the Earth's polar oblateness, or 2) use thrust to correct the injection errors. 

In this report only errors in semi-major axis and inclination will be considered. 
Also, all orbits will be assumed to be circular. 


ANALYSIS 

This section will deal with the analysis of the methods which could be used 
to obtain the sun synchronous precession rate. To reiterate, the two general 
methods are: 1) use a force normal to the plane of the orbit in such a manner 
that the desired value for the precession rate will be achieved, and 2) use 
thrust to correct either of the two injection errors that are being considered in 
this report (altitude and inclination). Also in this section, a relationship between 
altitude and inclination for a sun synchronous satellite will be obtained. 
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At the conclusion of this section a specific method for correcting the Tiros 
M precession rate will be given. The conclusion will be based on propulsion 
requirements only. 

Relationship Between Inclination and Altitude for a Sun Synchronous Satellite 


In this paragraph a formula which relates inclination and altitude for a sun 
synchronous satellite will be derived. 

First it is imperative to have a clear idea of what sun synchronous means. 
When the angle between the satellite's line of nodes and the projection into the 
equatorial plane of the line between the center of the Earth and the Sun has an 
average constant value, an orbit is said to be sun synchronous. This would 
mean that the angle C in Figure 1 would have an average constant value. 

In Appendix A it is shown that the average precession rate due to the 
Earth's oblateness is 




3 • (lA 

= "TJao^Vr/ 


cos I 


( 1 ) 


If it is assumed, as it will be throughout this report, that we have a near cir- 
cular orbit, we have 


T’q 


2 





and hence 


4>( deg/ day ) - -cr , J 20 -fl r e 2 r 7/2 cos I , 


where is a conversion factor. Putting in the constants gives 
0(deg/day) = - {2. 3817 x 10 13 ) r“ 7/2 cos I 

= +0.98S (the sun synchronous precession rate) 
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or 


I = cos -1 [-4. 136x 10‘ 14 r 7 2 ] (2) 

where r is given in nautical miles. Equation (2) now gives us a useful relation- 
ship between altitude and inclination for a sun synchronous satellite. Figures 2 
and 3 illustrate this relation in inclination versus altitude plots. Any point along 
these curves could be used for a reference orbit depending on the mission 
objectives. Note that atmospheric drag, and solar and lunar perturbations have 
been ignored here and will be throughout the report. 

Normal Thrust (Aid Oblateness) 


If a reference orbit is not obtained, two general ways to keep the precession 
at the correct value were given in the introduction. The first way mentioned 
was to use a thrust normal to the plane of the orbit in such a manner that it 
could change the precession rate. 

Our basic equation is 


2 M = H 


( 3 ) 


where, 

M - moment vector 

H = rate of change of angular momentum. 
In the case under consideration 


2 M = M p + M t 


where, 

— 

M p = precession moment due to Earth's oblateness 
M t = moment due to thrust. 
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for Sun Synchronous Rate 








It is assumed that the thrust comes solely from a force normal to the plane 
of the orbit and thus with the aid of Figure 4 can be written 


F n - F n | [-Sin I j + cos I k] 


- F n [-sin I j + cos Ik] 


(4) 


where, 


F n • 


Also, according to the nomenclature of Figure 4 


r 


r | [cos r) i + sin r) cos I j + sin 77 sin I k] 


[ A A A 1 

cos 17 i + sin 77 cos I j + sini7sinlkj 


( 5 ) 


where I r I - r , The moment due to thrust can now be determined from Equa- 
tions (4) and (5). 


M t - r x F N 


r F N {sin 77 i - cos 17 cos I j - cos 77 sin I k} 


( 6 ) 


The components of the precession moment that are determined in Appen- 
dix A are 


-rj f 1 - cos 2 t)] sin I cos I 

(7) 

= rj sin 2 rj sin I 

(8) 

Mpz = 0 

(9) 
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X (along line of nodes) 


Figure 4. Reference Coordinate System used in the Normal Thrust (Aid Oblateness) Paragraph 


where, 



r, 

& 
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The components of the moments in the Y and Z directions will be examined 
first. 


2 M v 


M py + 


M 


T Y 


- r J sin 2r) sin I - r F N cos 17 cos I 


H, 


( 11 ) 


where for reasons which will become apparent later, F N is so chosen that 





0 < 7} < n 
it < r) < 2rr 


( 12 ) 


If Equation (11) is now integrated over one orbit, the AH y resulting will be zero. 
A similar situation arises for the Z component. Thus, only the possibility of a 
change in H in the_X direction over an orbit is left. Since this is instantaneously 
perpendicular to H, the direction of H will be changed. 

Now the X component can be examined. Using Equation (3) again gives 


2 M x = Mp X + M 


TX 


= H 


(13) 


where 


H = ki x H 

= mr a 7j0sinl 


(14) 
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Rewriting Equation (13) using Equations (14), (6) and (7) gives 

mr 2 iyjj sin I = -rj [l - cos 277] sin I cos I + rF N sin 77 


Solving for 0, 


0 


J Cl - COS 2rj) COS I S * n ^ 

mr j) mr tj sin I 


(15) 


Assuming a circular orbit, i.e. 


fh = WT 

i - * ' 

Equation (15) can be integrated using Equation (12). The result is 

(t\ 2 4r 2 F„/m 

0( rad/orbit) = “^J 20 \^r) cos 1 * M sinl ’ 

or in degrees/day this becomes 


(16) 


(17) 


0( deg/ day ) 


1.9736 x 10* 
r */2 


0( rad/orbit) 


- +0.985° (for a sun synchronous satellite) 


( 18 ) 


Solving for the thrust to mass ratio 
M sin 

= J Z — 


V 


in 


I f 0.985 
\l.9736xlC 


3nJ 


20 


io* y? 
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or 



. . J4. 757 x 10” 1 1. 1502x 10 13 cosI 

sin I i — + 


(19) 


where r is in nautical miles. This equation indicates the amount of thrust that 
will be required at a given altitude and inclination to give a satellite of mass m a 
sun synchronous precession rate. 

F N /m versus I as determined from Equation (19) has been plotted for various 
altitudes in Figure 5. 

It should be noted that this is a continuous thrust that will last as long as the 
satellite is desired to remain sun synchronous. This implies that for missions 
of one to two years the total impulse required would be extremely large. 

Orbit Correction 

The second method mentioned to correct the precession, if the nominal orbit 
is not obtained, is to correct altitude or inclination to a sun synchronous value as 
given in Figure 2. 

Altitude Correction— To correct altitude, a low tangential thrust can be used 
to spiral in or out without significantly changing the eccentricity which is assumed 
to be initially near zero. In Reference 4, an equation is derived for total trans- 
fer time as a function of thrust and altitude change. It is 

yi j 

t t = ~ 1 ff L [l - ex P (-AV T /g I lp )] (20) 


where AV T can be calculated from 


AV t = )C/R^ - (21) 


Equations (20) and (21) were derived under the assumptions: 1) a nearly circular 
orbit, 2) a small thrust to mass ratio, and 3) a spherical Earth. 


It should be carefully noted that even though a spherical Earth was as- 
sumed, Equations (20) and (21) are still very good approximations for a non- 
spherical Earth. This can be shown by taking the gravitational forces for an 
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oblate earth given in Reference 5 and rotating them through the azimuth angle 
(azimuth angle here is defined from the equator) so that they will be in the same 
coordinates Reference 4 used. Two cases which can then be analyzed are: 1) 
when the ratio of Ar/r over an orbit is much less than one, and 2) when the ratio 
is not much less than one. When Ar/r << 1, the tangent gravitational force can 
be integrated over an orbit. This term integrates to zero, indicating that the 
Ar caused by this force is zero over an orbit since Ar is proportional to a tangen- 
tial force. If r is large enough so that Ar/r is not much less than one over an 
orbit, then the oblateness force is quite small in comparison to the central 
force term in the equations of motion and can thus be ignored. 

It should also be noted that by a method exactly identical to that performed 
in the "Normal Force" discussion, it can be shown that the only effect a tangential 
force has on H is that it changes its magnitude but does not add a precession 
torque. 

Using Equation (21) a plot of change in velocity versus change in altitude 
for various final altitudes has been made (see Figure 6). By using Figure 2, 
for a given injected inclination and altitude, a certain change in altitude can be 
found which will make a satellite sun synchronous with the injected inclination. 
This then can be used with Figure 6 to obtain the tangential change in velocity 
requirements. 

Inclination Correction— To analyze the correction of inclination errors, a 
method involving the linearization of the equations of motion will be used. Even 
though this will be done using the potential for a spherical Earth, the derivation 
for the inclination is still quite valid for a non-spherical Earth when we are con- 
sidering only a near circular orbit. This follows by considering the manner in 
which a plane change can be induced by a force. The required force is one that 
has a component normal to the orbit plane and (referring back to Figure 4) 
which points along H for -n/2 <tj< tt/2 and along -H for n/2 < 77 < 77/2 (or the 
reverse directions for a plane change in the opposite direction). If the force 
points along H in the interval 0 < 17 < 77 and along -H in the other interval, no net 
change in inclination occurs over an orbit due to cancellation. This means that 
for an oblate Earth which is symmetrical to the equatorial plane, no net change 
in inclination will occur due to the gravitational force since the force is of the 
type just mentioned. 

The equations of motion used in this derivation are: 


r = (r0 2 cos 2 <t> + r<£ 2 ) + F r /m - /V r 2 


( 22 ) 
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(23) 


5t (r 2 ^ cos 2 

5t (r 2 i) = 


*) 


(Fg/m) r cos <t> 


r ( V — r<9 2 sin 4> cos 0) 


(24) 


where r, 0 and 0are defined as in Figure 7. The reference orbit for th** lin- 
earization will be circular. The perturbed variables can be written 



Figure 7. Reference Coordinate System Used in the Inclination Correction Paragraph 
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0 


80 


(26) 


0 O + §0 0-80 




6 - n + h6 


6 = W 


(27) 


where r Q , 0 O and p refer to the circular reference orbit, and Sr, 80 and 8 9 are 
the perturbations with respect to this reference. Without any loss of generality 
let 0 O = 0. Equation (26) then becomes 


0 — 80 

0 - 80 

• • • • 

0 = 80 • (28) 


Substituting Equation (25), (27) and (28) into Equation (22) and neglecting terms 
second order or more gives 


Sr = r 0 n 2 + n 2 Sr + (r 0 + Sr) {2n && + ( 80) 2 } 

+ F r /m - n/r* + 2(^/r 0 2 ) 8r/r„ . (29) 


Since the reference orbit is circular, 


r 0 n 2 = M/r 0 2 • 


(30) 


If Equation (30) is substituted into Equation (29), and this result is divided by 
r 0 n 2 , we have 


1 8r 


n 


2 r, 





= F r /mr 0 n 2 (31) 
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In the exact same manner Equations (23) and (24) become 



By letting r - nt, Equations (31), (32) and (33) can be nondimensionalized. 


Also, all the non-linear terms will be dropped so that the system of equations 

can be solved. (The validity of this will be established later.) 
equations then becomes 

The system of 

d 2 /Sr\ Sr d 

dT ,(r 0 )-3 7;-2 dFW * * 

(34) 

d 2 d /8r\ Sr 

? (M)+2 4) tB ^ = B 

(35) 

d 2 Sr 

— (80) + S0 + C — = C 
a t o 

(36) 

where, 


A = F f / mr o n2 

(37) 

B = F e /mr 0 n 2 

(38) 

c = F V mr o ° 2 

(39) 
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The initial conditions and constraints on the system are: 

1) ?r/r 0 = B<p = $9 - 0 at r = 0 , and 

2) if A = B = 0, then 8r/r 0 = B9 = 0 . 

Also, since only low thrust to mass ratios, and radii less than 6,666 n.mi. (the 
maximum sun synchronous radius) are being considered, A, B and C « 1, and 
hence, terms second order or more in A, B and C will be ignored. 

Using these approximations and conditions, the solutions of Equations (34), 
(35) and (36) are 


$r/r 0 = 2Br 


(40) 


se 


2 t(A+ 3Bt) 


(41) 


80 - C {l - cos t + 2B( sin r - r)| 


(42) 


(It should be noted that if these solutions are substituted back into the original 
Equations, (31), (32) and (33), these equations will be satisfied if second or^er 
and above terms of A, B and C are neglected.) 

Since no perturbations in r and 6 are desired, let A = B = 0. This choice for 
A and B leaves only Equation (36) with a non-zero solution. Equation (36) becomes 


dl 

dr : 


— (80) + 80 = C . 


(43) 


If C is defined as 


C = 


J+ICI 

o < e < tt 

l-lcl 

77 < 9 < 2tt 


(44) 
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a 


Equation (43) will have a solution in 0 < \Q < -n like Equation (42), but in tt < 6 < 2n 
another solution arises since there will be new initial conditions (continuity at 
boundaries) and the sign of C has changed. If C is continued as a periodic func- 
tion of period 2n , the solution of Equation (43) will be 


80 



Vlcl 



(1 - cost) 
(1 + 3 cos t) 
(1-5 cos t) 


0 < 6 < 77 
TT < 6 < 27T 
2n < 6 < 3tt etc. 


The solution is represented by a plot in Figure 8. 

The relationship between inclination change and 80 is 


max | 80 1 = AI 


From Figure 8 and the relation just given 


AI 


mean 



(45) 


or 

^orbi*. = 4|C| 

= 4F7mr 0 n 2 

! 

AI orbit - 4F N /mr 0 n 2 (46) 

The substitution of F N ~ F^ is justified since 80 has been assumed to be small. » 

i 

* 


I 



£*0** K- , * 

w 


t 



Figure 8. Solution of Equation (43) 

The equation of prime concern here is Equation (45) rewritten as 


2£n_1_ 

mean “ IT m m T T 
2 

- ~v IT t t (47) 
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where it was assumed rn ~ iJI/t. This can also be rewritten as 


AV 


N 



( 48 ) 


where the substitution AV N = F N r T jm was used. To validate this substitution it 
will be assumed that the propellant mass is much less than the total initial mass 
of the satellite. Using Equation (48), Figure 9, which plots change in velocity 
versus mean change in inclination, was made. 

General Comparisons of Precession Correcting Methods 

In this paragraph the three methods just considered for the precession cor- 
rection will be compared. 

The method which uses a normal thrust to aid the Earth’s oblateness re- 
quires a large total impulse for missions of a half year or longer. This is due 
to the fact that a thrust is needed as long as a sun synchronous precession rate 
is desired. This will be illustrated for a Tiros M mission in the next paragraph. 

The two methods for orbit correction will now be compared. These two 
methods are: 1) use a tangential thrust to correct the altitude by spiraling in 
or out, and 2) use a normal thrust to cause an inclination change. The com- 
parison will be based on the total impulse required for each method. 

The two basic relationships that will be used are Equations (21) and (48) 
rewritten as 


and 



( 49 ) 




7TA I 


wean 



J£ 

-Ah 


( 50 ) 
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where cr 2 is the needed conversion factor from n.mi./sec. to ft./sec. since r 0 is 
given in nautical miles. Let C eff be defined as the ratio of the total impulse 
necessary to make a satellite sun synchronous using only an inclination change, 
to the total impulse necessary to make a satellite sun synchronous using only 
an altitude change, i.e. 

mAV N . 

Ceff = (assuming W p /W Q « l) 


or 


C*ff av t ’ (®*) 


If it is also assumed that Ah/r„ « 1, C eff reduces to 


” r 0 AI m.«n < de S> 

C eff ■ 57.3 Ah 


(52) 


where AI mean (deg)/Ah is the slope of the sun synchronous curve, Figure 3, if 
Ah and Al me<n are sufficiently small. 

To obtain a better understanding of C (f f and the variables which make it up, 
consider Figure 10. The solid curve represents the sun synchronous curve of 
Figure 3, and the dashed curve represents a constant precession rate error 
from the sun synchronous value. The point (r , , I , ) represents the injection 
point; (r 0 , 1 j ) represents the point on the sun synchronous curve which is 
reached by changing the altitude. Fro a Figure 10 it can be seen that a certain 
Ah and AI aeaii are associated with each injection point. 

Figure 11 shows C eff (which is in fact a penalty factor) as a function of 
altitude. This figure shows the interesting result that C eff is always greater 
than one. This means that the total impulse required to correct a given injection 
by correcting only inclination is always greater than that required for an altitude 
correction. 
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The conclusion that can be drawn from this is that it is always better, from 
a propulsion standpoint, to correct to a sun synchronous orbit with a tangential 
thrust (altitude correction). 

It should be pointed out, however, that the mission constraints may make 
this impossible. Referring again to Figure 10, where ( r, I ) represents the ref- 
erence orbit, one can see that an altitude error can arise if one only spirals in. 
In the next paragraph this will be illustrated further using parameters for 
Tiros M. 

Comparison of Precession Correcting Methods for a Tiros M Mission 


To illustrate what has been said concerning the various methods for cor- 
recting the precession rate error, the Tiros M mission will be considered. 
Tiros M is a near-polar, near-circular sun synchronous meteorological satel- 
lite. The nominal orbit for Tiros M is given in Table 1 (Reference 3). The 
1 cr injection error, given in terms of a precession rate error, is 0.029 degrees/ 
day (Reference 7). The total impulse required to correct the la injection error 
will be used to compare the three methods. 

So that the change in altitude and the change in inclination that are needed 
in Equations (28) and (41) can be found for a certain precession rate error, Fig- 
ure 12 and 13 have been made. Figure 12 plots change in precession rate in 
degrees/day versus change in inclination in degrees. Figure 13 plots change in 
precession rate in degrees/day versus change in altitude in nautical miles. 

The total impulse required to correct a precession rate error of 0.029 
degrees/day can be determined quite readily now. First it is noted that the la 
precession rate error could be caused by either an inclination error of 0.34 
degrees or an altitude error of 35 n.mi. (these are found from Figures 12 and 
13). Now using Figures 5 and 6 and Equation (48) the total impulse for each of 
the three methods can be obtained. The results are listed in Table 2. 

Table 2 illustrates quite clearly what was said in general earlier. The 
total impulse for correcting inclination is approximately 2.2 times as large as 
that needed to correct altitude; this is in agreement with the value given on the 

Table 1 

Tiros M (nominal orbit) 
inclination = 101.4° 
altitude = 750 n.mi. (approx.) 
weight of satellite = 670 lbs. 
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Table 2 


Total Impulse Required to Correct 
a Precession Rate Error of 0.029 deg. /day 


Normal Force (aid oblateness) 

138,700 lb.-sec./year 

Tangential Force (correct altitude) 

2,058 ib.-sec. 

Normal Force (correct inclination) 

4,571 lb. -sec. 


C eff graph. The total impulse needed to aid the precession due to obiateness 
for one year is approximately 67.4 times that needed to correct altitude. 

The above facts illustrate the basic conclusion: from a propulsion stand- 
point, it is best to use a tangential thrust, which corrects altitude, to get back 
to a point on the sun synchronous curve. 

As was mentioned earlier, one objection to using a tangential thrust is the 
error in altitude that might arise. Figure 14, which plots error in altitude (final 
sun synchronous altitude minus the reference altitude) versus the injected inclina- 
tion, illustrates this. For the la- injection error which corresponds to an in- 
clination error of 0.34 degrees, the altitude error would be approximately 35 
n.mi. An error of this magnitude could be tolerated by Tiros but for a satellite 
such as EROS, which has a very small allowable error in altitude, this would be 
far too much (Reference 6). 


EVALUATION OF PROPULSION SYSTEMS 

The purpose of this section is to compare propulsion systems that can be 
used to correct the Tiros M injection errors. Five propulsion systems are con- 
sidered for the mission. They are: 1) a 20 /xlb. ion engine system, 2) a 300 
Mlb. ion engine system, 3) two thermal storage resistojet systems, and 4) a 
2 lb. hydrazine system. Each system is assumed to be mounted so that the 
thrusters are tangent to the orbit (the configuration that was just shown in the 
preceding section to be the best from propulsion standpoint, see Figure 15). 

It should be carefully noted that each propulsion system contains two thrusters 
so that the satellite can be either spiraled in or out. Three aspects of each 
system will be compared. These are: 1) total weight of the system, 2) total 
time needed to correct the orbit, and 3) total power needed to operate the 
system. 
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Figure 14. Error in Altitude vs. Injected Inclination 






Orbit 



Figure 15. Thruster Configuration: Thrusters are Tangent to Orbit 


The propulsion system weight is made up of six parts: 1) the propellant 
weight, W p ; 2) the weight of thrusters, W th ; 3) the weight of feed systems, W fs ; 
4) the tankage weight, W t ; 5) the weight of the pov/er conditioner, W nc and 6) the 
telemetry weight, W tel . Written as an equation this becomes 


W 


ps 


w_ 


'th 


+ w fs + w, + w_„ + w 


pc 


'tel 


(49) 


However, 


W.h ♦ W 


fs 


pc 


+ w 


tel 


- ^0 = 


+ 


constant . 


(50) 




; 


and 


W 


t 


* K, W„ + K, 

1 D d 


(51) 


where Kj and K 2 are system constants (Reference 6). Putting Equations (50) 
and (51) into Equation (49) gives 


- (l + K,) W p + K 3 


(52) 


where K 3 - K 0 +K 2 . The propellant weight can be determined by 



(53) 


Since for all the propulsion systems being considered W p « W 0 , the loss of satel 
lite mass due to the expulsion of propellant is negligible. Thus total impulse 
can be written 


It = 



(54) 


Combining Equations (52), (53), and (54) gives 


W. 


p* 


Mi) 


V^T 

* X .p 


+ K, 


(55) 


Using Equation (55) and the constants for the various propulsion systems 
given in Table 3, propulsion system weight versus increment in velocity plots 
can be obtained (see Figure 16). Figure 16 shows that the 20 /xlb. ion engine 
system is the lightest after 28 fps change in velocity; while from 0 to 28 fps, 
the resistojet (I tp = 200 sec.) system is the lightest. The ranking of the systems 


z 

t. 
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Table 3 

Propulsion Systems (two thrusters) 


System #1 

20 Mlb. 

ion engine 


X *p 


5000 sec. 


K, 

= 

0 


k 3 

= 

16 lb. 


Power 


25 watts 

System #2 

CO 

o 

o 

cr 

ion engine 


X . P 

= 

5000 sec. 


Ki 

= 

0 


k 3 

= 

30 lb. 


Power 

= 

73 watts 

System #3 

thermal storage resistojet 


^ *P 

= 

200 sec. 


F 

= 

0.005 lb. 


*i 

= 

0.4167 


k 3 

= 

12 lb. 


Power 

- 

66 watts 

System #4 

thermal storage resistojet (cold gas) 



= 

100 sec. 


F 

= 

0.050 It. 


K, 

= 

0.4167 


k 3 

= 

12 lb. 


Power 

= 

negligible 

System #5 

2 lb. hydrazine engine 


*.p 


225 sec. 


Ki 

= 

0.26 



= 

16 lb. 


Power 

= 

negligible 
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Wp. (lb) 




Figuro 16 . Propulsion System Weight vs. Chong* in Volocity Required to Corroct Orbit 


34 






with respect to lightest average weight in the 0-100 fps AV range is: 1) 20 Mlb. 
ion engine system, 2) 200 sec., I sp resistojet system, 3) 2 lb. hydrazine system, 
4) 100 sec., I resistojet system, and 5) 300 /4b. ion engine system. 

The next comparison that will be made concerns the time required for a 
given correction. Since the satellite mass remains nearly constant, the total 
impulse can then be written 


It 



= Fr T 


or. 


T 


T 


w 0 av t 

gF * 


(56) 


Using the values given in Table 3 and Equation (56), Figure 17, which plots 
thrust time versus change in velocity, was constructed. Figure 17 clearly 
shows that the 2 lb. hydrazine system takes the least amount of time for any 
AV t required. The ion engine systems took considerably more time than any 
other system, e.g., the 20 /4b. ion engine system required approximately 486 
days for a 40 fps change in velocity and the 300 /4b. ion engine system re- 
quired approximately 32.2 days for the same change in velocity. Rating the 
systems according to least time required for a given change in velocity gives: 

1) 2 lb. hydrazine system, 2) 100 sec., I, p resistojet system, 3) 200 sec . I (p 
resistojet system, 4) 300 / lb. ion engine system, and 5) 20 /4b. ion engine 
system. 

The last comparison to be made is that of total power required. Using 
the powers listed in Table 3 (these were obtained from Reference 6) it can be 
seen that m terms of least power required the systems rate: 1) 2 lb. hydrazine 
system and 100 sec. I sp resistojet system (cold gas), which required almost 
no power, 3) 20 /4b. ion engine system, 25 watts, 4) 200 sec., I ^ resistojet 
system, 66 watts, and 5) 300 /4b, ion engine system, 73 watts. 

So that a comparison of all the propulsion systems can be more easily made, 
the Tiros M mission is again considered. Table 4 lists for each propulsion 
system the weight of the propulsion system, the thrust on time, and the power 
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Table 4 


Tiros M Propulsion System Requirements 

(Note: These are for an approximate precession rate error 
of 0.029 degrees/day) 


Propulsion Systems 


r T (days) 

Power (watts) 

20 /J-lb. ion engine 

16 

1208. 

25 

300 Ailo. ion engine 

30 

80.5 

73 

200 sec., I sp resistojet 

27 

4.8 

66 

100 sec., I resistojet (cold gas) 

42 

0.5 

negligible 

2 lb. hydrazine 

28 

0.012 
(1,043 sec.) 

negligible 


requirements needed to correct a precession rate error of 0.029 degrees/day 
(the 1 a error for Tiros M). A correction time of more than 30 days was con- 
sidered unacceptable for Tiros; thus, the ion engines are ruled out as possible 
systems for Tiros M. This leaves only the resistojets and the hydrazine system 
for further consideration. The final decision, on which system would be the best 
for Tiros M, rests on the attitude control system. It is possible that a thrust 
misalignment on the hydrazine system could cause disturbances in the attitude 
that could not be corrected by the control system. This problem, however, is 
not being treated in this paper. 


CONCLUSION 

Three methods for correcting a precession rate error, that can result from 
injection errors, have been considered in this report. They are: 1) use a force 
normal to the plane of the orbit to aid the oblateness force, 2) use a tangential 
force to spiral in or out to an altitude which gives a sun synchronous precession 
rate at the injected inclination and 3) use a force normal to the plane of the 
orbit to obtain an inclination which causes a sun synchronous precession rate at 
the injected altitude. The analysis showed that of the three methods considered 
tangential thrusting would be the best method from a propulsion standpoint. 

It should be kept in mind, however, that the final altitude may be different than 
the reference altitude if the injected inclination is different than the reference one. 
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Comparisons were made of the total system weights, the orbit correction 
times, and total power requirements of five propulsion systems (the systems 
include two thrusters to allow for either spiraling in or out). The five systems 
are: 1) a 20 /ulb. ion engine system, 2) a 300 /o.lb. ion engine system, 3) two 
resistojet systems, one with I = 200 sec. and thrust equal to 0.005 lb., and 
the other with I = 100 sec. and thrust equal to 0.050 lb., and 5) a 2 lb. hydrazine 
system. The comparisons of the five systems were based on the requirements 
ne' Jed to correct the lcr error for Tiros M (0.029 d^g./day). The ion engines 
required more time than was acceptable for the Tiros mission and were eliminated 
from further consideration. The weight of the 200 sec., I sp resistojet and the 
2 lb. hydrazine system were about the same. The hydrazine system (if fired 
continuously) required the least time for the correction. The cold gas resistojet 
and the hydrazine system required negligible amounts of power. An important 
factor that was not included in this paper was the effects of thrust misalignment. 
The attitude distrubances caused by a thrust misalignment must be compared 
with the control capability of Tiros before a final decision can be made. 
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APPENDIX 


PRECESSION DUE TO THE EARTH’S OBLATE NESS* 

The purpose of this appendix is to derive expressions for the moments and 
the orbit precession rate due to the equatorial bulge. 

From Reference 5 the forces due to the equatorial bulge are 

f r = J (3 sin 2 L - l) , 
f L = -^JsinLcosL , 

and 

Where 

J = 



Equation (Al) transforms tc components in the X, Y, Z coordinates of Fig- 
ure A-l as 


f x = j(5sin 2 L~l) cos \ cos L , 
f y = J (5 sin 2 L~ l) sin\ cosL , 


and 


f z - J (5 sin 2 L - 3) sinL 


(A2) 


•Derivation was performed by C. C. Barrett of the Auxiliary Propulsion Branch, Systems Analysis and Ion 
Propul nit..* Section* 
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and evaluating 


M r 


r x f 


gives 


M, 


px 


M 


P Y 



and 


M p2 = 0 . 


(A4) 


Noting that in terms of orbital variables co + £ and I 


— - cos (ai+ £) , 


— - sin (a: + £) cos I , 


and 


— = sin (&o + g) sin I' , 


(A5) 


Equations (A4) become: 

M px = rj [l - cos 2 (oj+ £)] sin I cos I , 
M py = rj sin 2(cv + £) sin I , 


and 


Mp Z 0 . 


4l 


(A6) 



All of M py and that part of M PX containing cos 2(o. + ~) have a net value of zero 
over one orbit, or as ('•- ± z) changes from 0 to 277. Therefore the net moment 
over an orbit acts about the X-axis and has a value of 


M pZ - -rj sin I cos I 





sin I cos I 


(A7) 


This moment produces a change in momentum given by 


H 


k i/7 x H 


i mr 2 £ \jj sin I 


(A8) 


Equating Equations (A7) and (A8) yields 


mr 2 ; 0 sin I 



Since this equation is valid for all inclinations, there follows 



(A9) 


For a nearly circular orbit £ ~ mA 3 and Equation (A9) can be rewritten as 


4 > 



cos I . 


(A10) 
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